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The dependence of pseudo-scalar masses and decay constants on the sea and valence quark masses is investigated 
in the pseudo-Goldstone boson sector of QCD with two light quark flavours. The sea quark masses are at present 
in the range |m s < m u d < \ m a whereas the valence quark masses satisfy |m S e« < rrivai < 2m sea . The values 
of the Gasser-Leutwyler low energy constants L4, L 5 , Lq and Lg are estimated. The computation is done with 
the Wilson-quark lattice action at gauge coupling j3 = 5.1 on 16 lattices. 0(a) effects are taken into account by 
applying chiral perturbation theory for the Wilson lattice action as proposed by Rupak and Shoresh. 



1. Introduction 

Monte Carlo simulations of QCD with dynam- 
ical quarks are done in most cases at relatively 
large quark masses (typically two quark flavours 
with m u d > m s /2). This makes the extrapola- 
tion to the physical point m u d ~ m s /20 rather 
uncertain. The extrapolation is done by using 
(PQ)ChPT - typically to NLO (1-loop) order. 
Estimates |1I2) show that one should perform sim- 
ulations in the range m u d < m s /4, m s /3 in or- 
der to see the expected logarithmic dependence. 
Matching the predicted functional dependence is 
a crucial test for Lattice QCD. 

Until recently the comparison between lattice 
data and ChPT was not satisfactory. In a recent 
paper we suggested that the agreement can 
be found when going to light enough dynamical 
quarks. Here we supplement some integration of 
the analysis of [3] and provide some further com- 
ments. 



2. Simulations 

For our simulation we used the algorithm de- 
scribed in 0] and further improved as in |5J and 
references therein. In our range of parameters we 
found that the cost for producing one indepen- 



dent gauge field configuration roughly goes as: 
Cost ~ F (am q y 2 n, 

where am q is the quark mass in lattice unit, f2 
is the number of lattice points and the factor 
F ~ 10 7 flop for the plaquette, but one order of 
magnitude smaller for and / w . 

We produced three sets of 0(1000) thermal- 
ized configurations for Nf = 2 unimproved Wil- 
son fermions at Vol= 16 , f3 = 5.1 and k 
0.1 7fj. 0.1765, 0.177. From these points we extrap- 
olated a value of ro/a at K cr = 0.1773(1) which 
is ro(K cr )/a — 2.65(7). This corresponds to an 
UV cutoff a = 0.189(5)fm ~ (1.04 GeV)" 1 , and 
to a physical volume L ~ 3fm. In those three 
points we found a ratio m q /m s equal to 0.68, 0.59 
and 0.35 respectively. Here m q represents the sea 
quark mass defined as M r — (m T ro) 2 , when the 
strange quark mass corresponds to M r — 3.1 [H]- 
For the lightest valence quarks the ratios m q al /m s 
become respectively 0.37, 0.34 and 0.17. Even 
with these small quark masses finite volume ef- 
fects are expected to be under control since we 
have always Lm^. al > 4.8. Of course this is payed 
with a very low UV cutoff, and one expects large 
lattice artifacts. These are taken into account in 
the analysis. 
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Figure 1 . RRf — 1 predictions against the data. 

3. Comparison with ChPT 

Our first goal is to confront the results of nu- 
merical simulations with the (PQ)ChPT formulas 
| till) . In order to cancel the Z-factors of mul- 
tiplicative renormalization, which in the case of 
a mass-independent renormalization scheme only 
depend on the gauge coupling and not on the 
quark mass, we considered ratios of quark masses 
(jn q A), pion masses (mab) and pion decay con- 
stants (/as)- Here A, B stand for the flavor in- 
dices of valence (V) or sea (S) quarks. 

If we assume that there are no lattice arti- 
facts, no NNLO corrections and we take for xs — 
2B m qS ^ ree ] eve j es timate, x| s * = Jr^p » 

then the ratios (here £ = — ) 

v s xs > 

£2 -.est 

RRf ee = 1 + ^f-aK - 1 - MOL 

JvvJss 32N s n< ! 

" (£+l) 2 m 2 vv m 2 s 16N s tt 2 

are non trivial and parameterless predictions of 
ChPT. This provides a strong check of how far we 
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Figure 2. RRn — 1 predictions against the data. 



are from the NLO ChPT regime. In figure ^ an d 
12 we superimpose the predicted functions to the 
data. The agreement is of course not complete 
(there are indeed 0(a) and NNLO effects), but 
the corrections are sub-dominant contributions. 
Moreover the agreement improves when the sea 
quark masses decrease. 

Encouraged by this results, we systematically 
compared our data with those ratios of pion 
masses and coupling constants which determine 
the Gasser-Leutwyler [7] coefficients L 4 , L 5 , L 6 
and Lg. Since we expected large lattice arti- 
facts in the data, the comparison was done with 
W(PQ)ChPT [B], including O(a) lattice artifacts 
in the effective continuum theory. Besides that, 
we also included the relevant contributions from 
NNLO ChPT 0. Although this involve many pa- 
rameters, one can obtain enough constraints from 
Partially Quenched simulations. The details of 
the fitting procedure and the results are described 
in [3]. Here we simply add that, following the 
analysis of ^HI> we can now include also 0(a 2 ) 
lattice artifacts in our fit of the pion mass ratios. 
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Figure 3. Comparison of the extrapolation of the Figure 4. Comparison of the extrapolation of the 

fit of RRf from /3 = 5.1 [3] with the original fit. fit of Rfvv from (3 — 5.1 3] with the original fit. 



It turns out that this does not add new parame- 
ters to the fit, but it amounts to a redefinition of 
the Wilson-ChPT coefficients W t . 

In order to have a consistency check of the sur- 
prisingly small lattice artifacts that we found, we 
compared the results of [3] with those obtained 
in an older simulation with larger lattice spacing 
(J3 = 4.68, a ~ 0.27fm) 11 . The comparison (fig. 
|3and0J shows quite small scale breaking. Larger 
scale breaking (~ 8%) are observed for the ra- 
tio m^/nip, which at fixed M r = 1.088 goes from 
0.613 (at = 5.1) to 0.664 (at /3 = 4.68). 

To summarize: we showed that it is possible 
to simulate, with reasonable costs, very light dy- 
namical quarks. Compensating O(a) effects by 
introducing 0(a) terms in the PQCh-Lagrangian 
is, in this case, a viable alternative to O(a) im- 
provement of the action. The observed 0(a) 
contributions in the light Goldstone boson sec- 
tor are surprisingly small, while NNLO are still 
important. The expected behavior predicted by 
PQChPT is already visible, although a quantita- 
tive determination of the LEC's still needs further 
simulations at smaller masses and lattice spacing. 
Most of the numerical calculations presented here 
have been done at the computers of NIC-Juelich 
and Zcuthcn. 
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